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We study a series of spectra measured in the superconducting state of optimally-doped Bi-2223 by 
scanning-tunneling spectroscopy. Each spectrum, as well as the average of spectra presenting the same gap, 
is fit using a microscopic model that takes into account the band structure, the BCS gap, and the interaction 
of electrons with the spin-resonance mode. After describing our measurements and the main characteristics 
of the strong-coupling model, we report the whole set of parameters determined from the fits, and we 
discuss trends as a function of the gap magnitude. We also simulate angle-resolved photoemission spectra 
using the parameters resulting from our fits, and compare with recent experimental data. 
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I. INTRODUCTION 

The two main single-electron spectroscopies, angle- 
resolved photoemission 1 (ARPES) and scanning tunneling 
microscopy 2 (STM), have been considerably improved dur- 
ing the last two decades, largely motivated by the need 
for ever more reliable data in the study of cuprate high- 
temperature superconductors. Thanks to a steady improve- 
ment in the sample quality, a rapidly growing body of high- 
quality spectroscopic data is now available for the cuprates, 
especially for compounds with easily-obtained clean sur- 
faces as the Bi-based family. These refined experiments on 
high-quality crystals seem now able to deliver the intrinsic 
line-shape of the one-electron spectra. Nonetheless, very 
few studies have so far undertaken a detailed line-shape 
study of the various spectral features, by means of a mi- 
croscopic model. In previous works, the STM data analysis 
remained largely qualitative, or based on phenomenological 
approaches. In a few cases, a BCS d-wave model includ- 
ing a realistic band structure, 3 and/or a phenomenological 
scattering rate, 4 turned out to be appropriate. Such cases 
are the exception rather than the rule: these simple mod- 
els will not capture, in particular, a structure ubiquitously 
present at energies above the superconducting gap, the "dip" 
feature. By analogy with phonon-related effects in classical 
superconductors, 5 ' 6 the dip, which has also recently been ob- 
served in pnictide superconductors, 7 ' 8 is generally attributed 
to a collective mode. 

The cuprates present, in their superconducting phase, 
a low-energy magnetic excitation known as the spin reso- 
nance. This excitation emerges below the superconduct- 
ing critical temperature, T c , as a strong enhancement of 
the spin susceptibility around the anti-ferromagnetic vec- 
tor Q = {n/a, Ti/a). First discovered in the inelastic neu- 
tron scattering response of YBa 2 Cu 3 6+x 9 ' 10 (Y-123) at 
an energy of fi 5 = 41 mey it has been observed since 
then in most cuprates, including the single-layered com- 
pounds HgBajCuO^g 11 and Tl 2 Ba 2 Cu0 6+5 , 12 the two-layer 
Bi 2 Sr 2 CaCu 2 8+5 (Bi-2212), 13-15 and the electron-doped 
material Pr 088 LaCe 12 Cu0 4 _ 5 . 16 The (71, n) resonance is 
detected at energies ranging from 10 to 60 mey roughly cor- 
related with T c following the dependence fi s 5.3k B T c . 17 
It is found in both underdoped and overdoped materials, 18 



and was also detected above T c in Y-123. 9 The resonance 
has been interpreted as the signature of a spin-1 exciton, a 
bound state below the continuum of electron-hole excita- 
tions, gapped by the superconducting pairing. 20 One of the 
theoretical questions, that remains unsettled so far, concerns 
the role played by spin fluctuations, and particularly by the 
(71, 71) resonance, in the pairing phenomenon. 21,22 

Being related to superconductivity or not, the (71, 7i) res- 
onance is a collective spin mode, which must somehow 
interact with the charge degrees of freedom, and induce 
renormalization and damping of the Bogoliubov quasipar- 
ticles in the superconducting state. Such effects certainly 
modify the one-electron spectra, and are in principle observ- 
able in single-electron spectroscopies. Although the strength 
of this interaction has been the matter of a controversy, 23 ' 24 
there is now enough evidence that peculiar signatures ob- 
served in photoemission, 25-28 tunneling, 29-33 and optical 
conductivity 34 ' 35 result from this interaction. Nevertheless, 
a firm consensus has not been reached yet: optical phonons 
are often present in the cuprates at nearly the same energies 
than the spin resonance, and discriminating unequivocally 
the effects of the two kinds of excitations has proven difficult. 
The spin resonance, being localized near (71, 71), naturally 
leads to an anisotropic scattering rate and a strong dip. 31 
But similar effects can also be induced by phonons, provided 
that the electron-phonon coupling is strongly anisotropic. 36 
A possible way of understanding the origin of the dip feature 
is to study its evolution as a function of doping, and to com- 
pare it with the doping-dependence of the spin resonance 
and of the phonons, both directly measured by neutron 
scattering. 

In a d-wave superconductor, one of the signatures of the 
coupling to a collective mode is a minimum, so-called dip, in 
the electron density of states (DOS), occurring at an energy 
E d , which is separated from the energy A p (of the coherence 
peak) by the mode energy, 37 E d = A p + Q. s . In contrast, fors- 
wave superconductors, the signature is a change of curvature 
of the DOS, leading to a peak in the DOS derivative. 38 ' 39 
Scanning tunneling microscopy/spectroscopy (STM/STS) 
probes the electron local DOS (LDOS) with sub-meV and 
atomic resolutions. 2 It enables to locally measure the d- 
wave gap in the electron excitation spectrum, as well as the 
dip, and to track their spatial variations in inhomogeneous 
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materials. 32 It is therefore an ideal tool to investigate the 
properties of the dip, and the relationship between the gap 
and the spin-mode energy. 

Among the cuprates, the three-layered compound 
Bi 2 Sr 2 Ca 2 Cu 3 O 10+5 (Bi-2223) is probably the best-suited 
material to perform such studies. Like other members of the 
Bi-family it can be easily cleaved, and offers atomically-flat 
surfaces for STM investigations. Bi-2223 has the highest 
optimal T c of the Bi-family, 111 K, a gap in the 30-60 meV 
range, and a very strong dip as revealed by tunneling 40 and 
photoemission 41-43 studies. 

In this paper, we study Bi-2223 STS spectra by means of 
least-square fits, with a microscopic model 37 containing the 
three ingredients required to reproduce the experimental 
data: the Van Hove singularity of a tight-binding band, a 
d-wave BCS gap, and a coupling to the (n,n) spin reso- 
nance. The results of similar previous studies were already 
reported. 31,32 Here, we describe our fits with full detail, pro- 
vide and discuss the complete set of parameters determined 
from the fits, and fit average as well as local individual spec- 
tra. Our motivation for performing direct fits to cuprate STS 
data is two-fold. The first aim is to demonstrate that, in spite 
of its simplicity, the model captures quantitatively the main 
characteristics of the tunneling spectra for optimally-doped 
Bi-2223: a V-shaped gap at low energies, tall coherence 
peaks and very pronounced dips, both features noticeably 
electron-hole asymmetric. Second, the good quality of these 
fits provides further evidence that the STM tunneling con- 
ductance measures the full electron LDOS, 3 rather than an 
effective quasiparticle DOS deprived of band-structure and 
self-energy effects. 44 

In Sec. II, we describe the growth and characterization 
methods of the samples used for this study, present the 
measurement method, and discuss the main features of the 
spectra. Section III is dedicated to the model used to fit 
the data. We use slightly different conventions than Ref. 37 
for defining the model parameters: for clarity we therefore 
describe and discuss our model in detail. We also explain 
our fitting method. In Sec. IV, we present our results and the 
trends in parameters deduced from the fits. We discuss the 
values of the most important parameters in Sec. V, compare 
with those obtained using other experimental techniques, 
and present simulations of ARPES intensities. Finally, we 
summarize the main results and the implications of the 
present study in Sec. VI. 



II. EXPERIMENTAL DETAILS 

A. STM measurements, sample growth and characterization 

The LDOS measurements were performed by means of a 
home-built STM with ultra-high vacuum environment, and 
3 He base temperature. 45 Electrochemically-etched iridium 
tips served as the ground electrode, and a bias-voltage V 
was applied to the sample. Therefore, negative (positive) 
bias refers to occupied (empty) sample states. Differential- 
conductance spectra, dl/dV, were directly acquired using a 
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FIG. 1. Exponential decay of the tunnel current with increasing 
the relative tip-sample distance z, for one of our typical junctions. 
These data were obtained in a measurement with the current 
feedback-loop open. Insert: same data on a logarithmic scale. 



lock-in technique. 2 All the dl/dV measurements reported 
here were performed at fixed tip-sample distance, deter- 
mined by regulation current and voltage of 0.6 nA and 
600 mV, and a lock-in excitation amplitude of 2 mV 

The samples were cleaved at room temperature at 1- 
5 x 10~ 9 mbar pressure, and cooled down to 2 K in 10 
hours time. In this way, high-quality tunnel junctions were 
typically obtained, as revealed by the data shown in Fig. 1 . 
In a high-quality tunnel junction, the current depends ex- 
ponentially on the relative tip-sample distance, z, namely 
I oc e~ 2KZ . The decay constant k is related to the apparent 
barrier height </> = H 2 k 2 /(2m). Our tunnel junctions have 
typically <f> = 3-4 eV, much larger than the regulation volt- 
age. We started all our runs of measurements by performing 
tests like the one shown in Fig. 1. 

The Bi-2223 crystals were grown by the traveling-solvent 
floating-zone method, as described in Ref. 46. In order to 
achieve optimal doping, the crystals were annealed dur- 
ing 10 days at 500°C under an oxygen partial pressure of 
20 bar. This thermal treatment produced a sharpening of 
the superconducting transition with respect to the as-grown 
condition. Figure 2(a) shows the low-field susceptibility ver- 
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FIG. 2. Magnetic properties of an optimally-doped Bi-2223 crystal 
from the same batch as the one used for STM studies, (a) Real part 
of the magnetic susceptibility as a function of temperature, reveal- 
ing a T c = (110.5 ± 0.5) K. Measurements were performed with an 
ac field of 0.1 Oe and 970 Hz. (b) Temperature evolution of the 
magnetization in field-cooling processes at various applied fields. 
The magnetization is shown normalized by the low-temperature 
saturation value M n . 
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FIG. 3. The complete experimental dataset considered in this study. The raw Bi-2223 tunneling conductance spectra were grouped 
according to their peak-to-peak gap A , and were normalized in order to have the same spectral weight in the range of the figure (black 
curves). The red curves are the average spectra, and the shaded blue curves show the standard deviation of the distribution of tunneling 
conductances at each bias. The vertical bars indicate a weak feature, possibly related to the Van Hove singularity of the inner-layer band 
(see Sec. III). 



sus temperature for a typical optimally-doped (OPT) sample, 
with a critical temperature T C = (110.5±0.5)K. 47 We con- 
sidered the peak position in the temperature derivative of 
the low-field susceptibility and magnetization data as the 
criteria to determine T c . Both determinations yielded the 
same T c . The transition width, estimated as the FWHM of 
the susceptibility and magnetization peaks, typically ranges 
between 0.6 and 2 K from sample to sample. 

The structural and superconducting properties of OPT 
crystals of the same batch than the ones studied here are 
reported in Refs 47-50. X-ray diffraction measurements 
have revealed the purity and the high crystalline order of the 
samples. 49 Electrical resistivity data also revealed a single 
and sharp superconducting transition. The sample-growing 
and thermal-treatment parameters were optimized, in order 
to suppress the intergrowth of the Bi-2212 phase in the 
three-layered compound. For the case of the OPT samples 
studied here, Bi-2212 intergrowth, if present, represent less 
than 1% of the sample volume. 47 



B. Data statistics and systematics of spectral features 

As largely documented in the literature, Bi-based cuprate 
samples, even with a sharp superconducting transition, 
present inhomogeneous spectroscopic properties on the sur- 
face at the nanoscale. 2 In a previous work, 32 we mapped 
the local dl/dV curves in order to study the nanoscale vari- 
ations of typical spectral features. In particular, we found 
that the local value of the superconducting gap measured by 
A p presents spatial variations in register with the crystalline 



structure. 

The present paper is focused on studying relevant spectral 
features as a function of A p . Since results in several OPT 
samples are similar, for the present study we considered 
150 spectra measured at different locations of a given OPT 
sample. The local dl/dV curves were sorted by half the 
peak-to-peak energy separation, A p = (A + — A~)/2, where 

(A~) is the energy of the coherence peak at positive 
(negative) sample bias. In order to calculate the average 
spectra for a given A p , individual curves were first nor- 
malized, such that all have the same integral in the range 
[—300, +300] mV Figure 3 shows the set of individual and 
average spectra for A p between 36 and 54 meV Also shown 
is the standard deviation to the average, as a function of 
bias. For all A p , the standard deviation is of the order of 
5%, except close to the coherence peaks where it increases 
up to 20%. 

The prominent spectral features of Bi-2223 are evident 
from the average differential conductance curves shown in 
Fig. 3. All spectra present a d-wave like spectral shape at 
low bias, very well developed coherence peaks at the gap 
edges, and a dip-hump structure at energies larger than A p . 
These features, as we discuss in the following paragraphs, 
systematically evolve with A p . Another remarkable property 
of the Bi-2223 spectra is a strong electron-hole asymmetry, 
characterized by a greater spectral weight for the occupied 
sample states. This is manifested by higher coherence peaks, 
and an enhanced conductance background at negative bias. 
Detected in most high-T c 's, this phenomenology contrasts 
with the symmetric spectra of classical superconductors. The 
line-shape of the coherence peaks also follows a monotonic 
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trend with A„, when considering its evolution for either 
negative or positive bias. The peaks sharpen and become 
taller on decreasing their energy location with respect to 
the Fermi level. 

The electron-hole asymmetry monotonically evolves with 
A p . Firstly, the conductance background for the occupied 
states becomes steeper on decreasing the gap. The trend is 
particularly evident when considering data in a bias range 
one order of magnitude larger than A p . This evolution 
may be due to a variation of the correlation effects with the 
pairing strength. 44 Since our microscopic model does not 
consider such correlations, we will restrict our fits to the 
bias interval [—150,-1-150] mV Secondly, the asymmetry in 
the height of coherence peaks systematically increases when 
decreasing A„. We will discuss the role of the Van Hove 
singularity in this phenomenology in Sees. Ill C and IV. 

The evolution of the dip is of particular importance for 
this work. This feature consists in a local minimum of 
conductance at E d , in the vicinity of the coherence peaks, 
outside the gapped region. In the case of Bi-2223, the dip 
feature is strongly developed in comparison to other Bi- 
based cuprates. The line-shape of this feature is noticeably 
electron-hole asymmetric, the dip being stronger in the 
occupied part of the spectrum. As we will discuss in the 
following sections, another important trend is the evolution 
of the energy separation E d — A p = Q. s with A p : fi s decreases 
on increasing the gap. 

Finally, the low-energy conductance is similar for all val- 
ues of A p . Close to the Fermi level, the spectra are electron- 
hole symmetric, and present a slightly rounded shape. A 
strictly V-shaped low bias conductance is expected for a d- 
wave superconductor at zero temperature. However, some 
rounding off and a finite zero-bias conductance are expected 
when considering thermal and measurement broadening, 
the latter due to the finite amplitude of the ac lock-in exci- 
tation, or to electronic noise. Our fitting procedure takes 
these two experimental factors into account. 



III. MICROSCOPIC MODEL AND FITTING PROCEDURE 

A. STM tunneling conductance and LDOS 

The theory of tunneling in superconductors 38 was histori- 
cally meant for planar junctions involving classical s-wave 
superconductors, with a structureless normal-state DOS. 
On the other hand, the theory of Tersoff and Hamann 51 
for the STM was not developed for superconductors. The 
approach of Ref. 51 is nevertheless easily extended to de- 
scribe STM tunneling into superconductors, 2 and leads to 
the well-known paradigm that the differential conductance 
is a measure of the thermally-broadened electron LDOS: 
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N(r, co) is the sample LDOS at the position, r, where the 
STM tip is located, /' is the derivative of the Fermi function, 



and M is a tip-dependent constant. We have considered 
in Eq. (1) an additional broadening, in order to take into 
account the finite experimental resolution. This is modeled 
by a Gaussian function g a , with half width at half maximum 
a. There are several sources of extrinsic broadening in 
the experiment, among which the most obvious are the 
lock-in ac modulation, and the averaging of several similar 
spectra (see Sec. II). The strict proportionality of dl/dV and 
N(r , eV) is achieved in the limits of zero temperature and 
<j = 0, where both — /' and g a become delta functions. 

When the work function of the junction is large com- 
pared with the typical energies of interest (in our case, 3 eV 
compared with 0.1-0.2 eV see Fig. 1), it is reasonable to 
assume that all Bloch waves decay exponentially outside 
the sample surface, with the same energy-independent de- 
cay constant k. The LDOS at the tip location r reduces 
then to N(r,co) oc exp(-2K-z)JV(o>), with N(co) the bulk 
two-dimensional DOS, assumed translation-invariant in the 
(x,y) plane. Consistently, the current must decrease expo- 
nentially with z, as confirmed in Fig. 1. The z-dependent 
factors, which are irrelevant in spectroscopic measurements, 
are lumped into the constant M of Eq. (1), where N(r, co) 
can thus be replaced by JV(<u). We then calculate the two- 
dimensional DOS as the Brillouin-zone integral of the elec- 
tron spectral function: 



AT(cu)=-^](-l/7r)ImG 11 (fe ) w). 



(2) 



Here, G n is the first component of the one-electron Green's 
function in the Nambu particle-hole representation, and JV 
is the number of k points in the two-dimensional Brillouin 
zone. In a superconductor characterized by a BCS gap A fc , 
and inelastic scattering processes, it is convenient to write 
G n in terms of the 2x2 matrix self-energy co), in the 
form 
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(3) 



S 11 (k, co) and £ 22 (k, co) describe the renormalization and 
damping of the Bogoliubov quasiparticles in the electron 
and hole branches, respectively, while the "anomalous" self- 
energy £ 12 (k, co) describes inelastic-scattering effects in the 
pairing channel. 52 The explicit form that we use for S is 
provided in the next subsection. If £ = 0, Eq. (3) reduces 
to the classical BCS Green's function, with E, k = e k — \x the 
noninteracting electron dispersion, \i the chemical potential, 
and r a phenomenological constant scattering rate. 53 We 
use a tight-binding model for the band e k , involving up 
to five-neighbor hopping, which reads (setting the lattice 
parameter a = 1): 



e k = ^t(|r|)e" £r = 2t 1 (cosfc x + cosfc y )+ 

r 

4t 2 cos k x cos ky + 2t 3 (cos 2k x + cos 2k y )+ 
4t 4 (cos 2k x cos k y + cos k x cos 2k y )+ 

4t 5 cos2fc x cos2fc y . (4) 
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Note that our conventions for the signs and magnitudes of 
the hopping amplitudes t, differ from those used in Ref. 37. 
For the BCS gap, we assume a pure d x 2_ y i symmetry, i.e., 
A fc = A (cosfc x — cosfc y )/2. 

Bi-2223 has three Cu0 2 layers per unit cell, and therefore 
displays three bands at the Fermi level. 54 We nevertheless 
use a one-band model. This model turns out to be appro- 
priate for fitting the spectra in the range [—150, +150] mV, 
and has the advantage of reducing the number of adjustable 
parameters. Recent photoemission studies suggested that 
the Bi-2223 band triplet is actually divided into a doublet of 
nearly degenerate bands, attributed to the two outer Cu0 2 
layers, and a single band at lower energy, associated with 
the inner Cu0 2 layer. 55 ' 56 The inner-layer band is seen be- 
low the doublet, at ~ 70 meV in the nodal direction, and 
at an extrapolated energy of ~ 200 meV in the anti-nodal 
direction. Therefore, a one-band approximation focusing on 
the doublet is sufficient to fit our spectra, since the latter 
carries 2/3 of the spectral weight, and the Van Hove singu- 
larity of the inner-layer band lies ~ 200 meV below that of 
the doublet, just at the border of our measurement window. 
At low energies, we therefore expect that the modifications 
induced in the theoretical STM spectra by using a multi- 
band description would be marginal. That said, we point 
out that the average spectra displayed in Fig. 3 systemati- 
cally present a weak structure between —215 and —225 my 
which might be the signature of the inner-layer Van Hove 
singularity. This is supported by the fact that the structure is 
only seen at negative bias. Although a definitive assessment 
is not possible at this stage, this observation confirms that 
possible multi-band effects are likely to be small. 



B. Bogoliubov quasiparticles coupled to spin fluctuations 

The theoretical investigation of the coupling between Bo- 
goliubov quasiparticles and spin fluctuations began with 
the study of the superfluid transition in 3 He (see Ref. 57 
and references therein), and was revived with the discov- 
ery of high- T c superconductors. 58-50 The minimal model to 
describe the effects of this coupling is 



the spin response is dominated by a resonance at energy fi s , 
near the antiferromagnetic vector Q = (71, n): 



Xs (q,m n ) = W s F{q) 



de 



if2„ - e' 



(6) 



We choose the real functions F(q) and /(e) such that 

F W = J™ deI ( £ ^> = L Thus > w s stands for the 
momentum and frequency integrated spectral weight of the 
resonance: 
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The function F(q) is Lorentzian-like, peaked at q = Q, with 
a half width at half maximum Aq : 



F(q) 



F 



sin 2 (2^) +sin 2 (^) + (Aq/4) 2 ' 



(8) 



The constant F ensures the normalization of F(q). In 
Ref. 37, the resonance was assumed to be sharp in energy, so 
that its energy distribution was /(e) = 5(e — fij — 5(e + fi s ). 
Indeed, neutron scattering measurements suggest that the 
resonance is resolution-limited in Y-123. 9 We have used a 
slightly more general form, 



/(e) = /„ [L r _(e - n s ) - L Fs (e + n,)] , 



(9) 



where L r (e) = (r/7i)/(e 2 + T 2 ) is a Lorentzian, and / 
ensures the normalization of /(e). The form (9) allows us 
to take into account a finite lifetime t 5 ~ T" 1 of the spin 
excitation. Neutron scattering experiments indicate that the 
resonance is somewhat broader in Bi-2212 (Ref. 13) and 
Bi-2223 (Ref. 63) than in Y-123, and would be consistent 
with T s 4—8 meV Alternatively, Eq. (9) may be regarded 
as a crude way to account for the experimentally observed 
dispersion of the resonance, 64 which would broaden the 
sharp mode into a band of width F s . 

It is convenient to write the Matsubara Green's function 
in Eq. (5) using the spectral representation: 



E(fc ) co) = --^]-2]^ s (q ) in„) 

1 1 ;n„ 

x %{k-q,ia> n -m n )\ iijn ^ +iQ+ . (5) 

% is the 2x2 Nambu-BCS -Matsubara Green's function in 
the absence of coupling, x s IS trie spin susceptibility, ico n and 
ifi n are the fermionic and bosonic Matsubara frequencies, 
respectively, f3 = (fc B T) _1 is the inverse temperature, and 
g is the coupling parameter. Equation (5) can be obtained 
from perturbation theory in the electron-spin coupling; 61 
it can also be viewed as a simplified, non-self-consistent, 
version of a conserving strong-coupling theory. 62 

Following Ref. 37, we use a separable phenomenological 
expression for % s in the superconducting state. In the energy 
range of interest, namely below ~ 150 meV, we assume that 



de 



i<d„ — e 



(10) 



Taking into account the phenomenological scattering rate T 
appearing in Eq. (3), the spectral function can be expressed 
in terms of Lorentzian functions as 



Mk, e) = u k L T {e -E k ) + v fc L r (e + E k ), (11) 



where 



Uk= 2 



At ^ej, v fc = ii-u fc , (12) 

E k J 



V E k 



and E k = \J £, \ + A 2 .. Inserting Eqs. (6) and (10) into 
Eq. (5), one can perform the summation over Matsubara 
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frequencies using standard techniques, 65 and using Eqs. (9) 
and (11), one obtains the self-energy on the real-frequency 
axis: 

a 2 

t{k,(o)= — F(q )[u fc _ tJ B(to, £ fc _ q ) 
i 

+ v k _ q B(co,-E k _ q )]. (13) 

We have introduced the dimensionless coupling constant 
a 2 = (g/A) 2 W s I , with A a characteristic energy scale of the 
model, which we take as the nearest-neighbor hopping t 1 . 
For definiteness, the function B(co,E) is derived in the Ap- 
pendix. The self-energy (13) is a convolution in momentum 
space, and can therefore be efficiently evaluated numerically 
on dense k -point meshes using fast Fourier transforms. 

The use of a dimensionless coupling constant a incorpo- 
rating the spectral weight of the resonance, instead of the 
original coupling g like in Ref. 37, is more convenient for 
our purposes. The strength of the self-energy effects on 
low-energy Bogoliubov quasiparticles — hence the strength 
of the dip in the tunneling spectrum — is governed by the 
product g 2 W s , and W s strongly depends on the momentum 
width of the resonance, Aq. In the original formulation, 
both g and Aq strongly affect the dip strength, whereas, in 
our formulation, the dip strength is controlled mostly by a, 
and weakly depends on Aq. Since the model parameters 
will be determined by least-square fits to the experimental 
data, one may expect a simpler landscape if one avoids 
that different parameters have the same influence on the 
theoretical spectrum. 

C. Discussion of the model 

The model has fifteen parameters, including the multi- 
plicative constant M in Eq. (1), thirteen of which are deter- 
mined by least-square fitting. The two fixed parameters are 
the experimental temperature T = 2 K, and the Gaussian 
broadening a in Eq. (1), set to 4 meV The latter value rep- 
resents the broadening effect of the lock-in modulation, and 
the average of several spectra to reduce the noise and pro- 
duce representative spectra, as shown in Fig. 3. Additional 
broadening is provided by the Dynes parameter T, which 
accounts phenomenologically for impurity scattering, and 
varies between 1 and 9 meV in our fits. The main parameter 
is A , which sets the scale of the d-wave gap. In this model, 
the d-wave gap is not assumed to result self-consistently 
from the coupling to spin fluctuations, in the fashion of the 
Eliashberg theory for phonons. Instead, it is assumed that a 
d-wave gap is preexisting, and that the coupling to the spin 
resonance comes as a perturbation. The effect of the spin 
resonance on pairing is discussed further below. 

The band parameters t x _ 5 and the chemical potential /x 
determine the noninteracting DOS, the band filling, and the 
noninteracting Fermi surface. The property of the band on 
which the tunneling spectrum is most sensitive is the energy 
of the Van Hove singularity (VHS), equal to the dispersion at 
the M point (7i,0), i.e., £ M = 4(-t 2 + 1 3 + t 5 )-/i. The effect 
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FIG. 4. Evolution of the theoretical tunneling spectrum with vary- 
ing five relevant parameters of the model. The base parameters are 
(all energies are in meV) M = 1, r = 2, t H = (-250, 80, 0, 0, 0), 
H = -320, A = 40, n s = 40, T s = 5, A, = 1.5/a, a = 0.6, 
and they correspond to the middle red curve in each series . The 
k -point mesh contains 1024 x 1024 points. In (a), the chemical 
potential is varied to move the Van Hove singularity with respect 
to zero energy from —80 to +80 meV; the thin lines show the cor- 
responding noninteracting DOS. (b), (c), and (d) show the result 
of varying the mode energy f2 s , the mode q -space width Aq, and 
the coupling a, respectively. The effect of changing r s is illustrated 
by the thin curves in (b). 



of varying £ M is illustrated in Fig. 4(a). Four main trends 
can be observed: the position of the VHS affects (1) the 
difference in the heights of the two coherence peaks, (2) the 
electron-hole asymmetry of the dip, (3) the overall height 
of the coherence peaks, and (4) the half peak-to-peak gap 
A p . (1) results from the opening of the superconducting 
gap, which pushes the VHS further down if initially located 
at negative energy, and further up in the opposite case. 
The coherence peak closest to the singularity thus carries 
more spectral weight. One notices that some electron-hole 
asymmetry remains when £ M = 0, even if the weight of 
the VHS is evenly distributed among the two coherence 
peaks. The reason is that the noninteracting DOS breaks 
electron-hole symmetry, since a t 2 # produces an excess 
weight at negative energies. (2) is a consequence of the dip 
being reinforced by the VHS, and therefore being stronger 
at negative energy if £ M < 0, and vice-versa. 31 (3) reflects 
the overall contribution of the VHS to the coherence peaks, 
so that the latter are tallest when £ M = 0. Lastly, (4) is a 
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consequence of the d-wave symmetry of the gap, implying 
that the maximum gap on the Fermi surface is generally 
smaller than A , and decreases as the Fermi surface moves 
away from the M point. 

The remaining four parameters characterize the spin reso- 
nance. Changing the resonance energy fi s has mainly three 
effects, as shown in Fig. 4(b): (1) the dip minimum moves 
with respect to the closest coherence peak, (2) the height 
of the coherence peaks changes, and (3) the gap size varies. 
(1) indicates that the damping of Bogoliubov excitations 
is strongest near the energy A + fi s : excitations near this 
energy can easily decay by emitting a (ti,k) spin mode, 
since there are many final states available near the gap-edge 
energy A . In order to figure out the typical energy depen- 
dence of the electron scattering rate, — ImE 11 (fc, co), one 
can use the simple limit of Eq. (A8) for the B function, and 
average Eq. (13) over the Brillouin zone, to arrive at 

- ImE n (w) = |(«ti) 2 [0(w - n s )N (co - n,) 

+ 6(-a>-n s )N (co + n s )], (14) 

where IV (a>) is the d-wave BCS DOS. Since JV (a>) has 
coherence peaks at \co\ w A , -ImS n (w) peaks at co fas 
±(A + Q. s ), leading to the occurrence of dips in the full 
DOS at these energies. The actual scattering rate has a 
marked k -dependence, 37,39 which is washed out in Eq. (14). 
The momentum dependence is due to the resonance being 
confined around (71, tz), and approximately amounts to in- 
serting a k -dependent amplitude in the right-hand side of 
Eq. (14). The latter still captures the correct qualitative en- 
ergy dependence. The effect (2) — increase of the coherence 
peaks height with increasing Q s — can also be understood 
with the help of Eq. (14): the scattering rate vanishes in 
the energy range \co\ < Q s , where Bogoliubov quasiparti- 
cles cannot decay, and thus the coherence peaks are not 
broadened if £l s > A . The broadening of the coherence 
peaks is maximum when the peak in — ImS n occurs at the 
energy A p , i.e., when fi s — » 0. The trend (3) has its origin in 
the renormalization of the quasiparticle energies, encoded 
in the real part of the self-energy. At low energy, we have 
ReE n (o)) » —Xco, so that the energy levels are renormal- 
ized by a factor 1/(1 + A). Performing the Brillouin-zone 
average as above leads to 



A = (at 1 ) 2 



de 



(15) 



If the anti-nodal Fermi crossing is not too far from the M 
point, the scale A p of the gap in the theoretical spectrum 
is given to a good approximation by A (l + * M )/(1 + A M ). 
A M is the renormalization at the M point, that is typically 
30% larger than the Brillouin-zone average A, and *I/ M = 
ReS 12 (M, 0)/A gives the contribution of the anomalous 
self-energy to pairing. ^ M is a decreasing function of fi s ; 
this can be seen by estimating the Brillouin-zone average of 
ReS 12 (fc,0)/A fc , which gives 



* » (at a ) 2 



de 



s{e + n s y 



Equation (16) is accurate in the limit Aq — » 0, and for a 
electron-hole symmetric band. From the Q s -dependence of 
* alone, one would expect a decrease of A p with increasing 
fi s . However, this is over-compensated by the faster decrease 
of A with fi s [see Eq. (15)], so that the net result is a slight 
increase of A p with increasing fi s , as seen in Fig. 4(b). 
Figure 4(b) also illustrates the effect of increasing the energy 
width F s of the spin resonance: that is a broadening mostly 
confined to the neighborhood of the dip minimum. 

Figure 4(c) shows that the momentum width Aq of the 
spin resonance mainly controls the gap renormalization. If 
Aq < a -1 , * M wins over A M , and consequently A p > A , 
while the opposite occurs if Aq > a -1 . The precise value of 
Aq for which this change of behavior takes place depends 
on the other model parameters. In the limit Aq — » 0, the 
renormalization parameters at M can be evaluated as 



[at,) 2 



(Aq = 0), (17) 



and clearly ^ M > A M in this case (£ M = y^M + Aq). In 
the opposite limit, Aq — » 00 or F(_q) = 1, the self-energy 
becomes momentum independent. We then simply get 



A M — A, 



*M = 



(Aq = 00), 



(18) 



(16) 



such that the gap is reduced by a factor 1/(1 + A) in this 
case. The vanishing of * M is a consequence of the d-wave 
symmetry of A fc . It is important to remember that the varia- 
tion in Fig. 4(c) is performed at fixed spectral weight of the 
resonance. In other words, the bare coupling g decreases as 
Aq increases, so as to keep the a coupling unchanged. 

Finally, changing the dimensionless coupling a produces 
a three-fold effect, as seen in Fig. 4(d): increasing a (1) 
digs the dip, (2) reduces the gap, and (3) lowers the coher- 
ence peaks without broadening them. (1) reflects the fact 
that the scattering rate is proportional to a 2 [Eq. (14)]. (2) 
results because Aq > a -1 in Fig. 4(d), so that * M < A M , 
as discussed previously; since both * M and A M are propor- 
tional to a 2 , the quantity (1 + * M )/(1 + A M ) decreases with 
increasing a. Finally, the effect (3) is due to the fact that cou- 
pling to the spin resonance removes spectral weight at low 
energies, and transfers it to higher energies, not only over 
intermediate energy scales corresponding to the "hump", 
but over a larger energy scale, up to the band edges. This is 
demonstrated in Fig. 5, which shows the data of Fig. 4(d) on 
an expanded energy scale. The transfer of spectral weight 
from low to high energies explains why the coherence peaks 
are lowered with increasing a. We further note that part 
of the spectral weight removed in the dip is also pushed to 
lower energies and, under certain conditions, can lead to 
shoulders on the sides of the coherence peaks. 

We close this section with a few general remarks. Unlike 
in the strong-coupling theory for conventional superconduc- 
tors, the model that we are using does not allow us to find a 
simple relationship between its parameters and the location 
of prominent features in the theoretical tunneling spectrum. 
The rule is that each feature is controlled in different ways 
by several parameters. For instance, the peak-to-peak gap 
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FIG. 5. Same data as Fig. 4(d) on an expanded energy scale 
covering the whole bandwidth. Part of the low-energy spectral 
weight is transferred to the band edges. 



A p depends on A , £ M , fi 5 , Aq, and a, as illustrated in 
Fig. 4. A direct readout of the parameters from inspection of 
the spectra is therefore possible only qualitatively, and fitting 
is required for extracting accurate values. Yet, there is one 
peculiarity of the spectrum, which can be attributed to a sin- 
gle source: the electron-hole asymmetry — of the coherence 
peaks height, dip strength, and conductance background — is 
entirely due, in this model, to the electron-hole asymmetry 
of the underlying normal-state DOS. Therefore, the qualita- 
tive inference of an asymmetric normal DOS can be made by 
direct inspection of the spectra. Finally, we emphasize that, 
although the contribution of the anomalous self-energy to 
pairing is always positive, the experimental magnitude A p 
may eventually turn out to be smaller than A , due to the 
normal self-energy renormalization. Namely, the coupling 
to the spin resonance contributes positively to pairing, but 
can nevertheless reduce the preexisting gap. 



In a first step, we have focused our attention to the av- 
erage spectra displayed in Fig. 3. For each of them, we 
searched a set of parameters which (i) is a minimum of the 
least-square function, 67 (ii) corresponds to a noninteracting 
band with plausible hopping parameters, and the properties 
shared by all Bi-based cuprates — band minimum at the F 
point, band maximum at the (n, n) point, VHS at the M 
point — and (iii) provides a hole-like Fermi surface centered 
at (7i, 7i). We have not constrained the band parameters in 
order to fix the carrier density, because the local density cor- 
responding to each spectrum is unknown, and the nominal 
density deduced from the bulk T c is not very meaningful in 
heterogeneous samples. The actual procedure has been to 
search good parameters for the spectra with A p = 36 and 
54 meV (this required a bit of trial and error), and then to 
use interpolations between these parameters as seeds to fit 
the average spectra with intermediate gaps. The fits were 
restricted to the energy window \co\ < 150 meV 

In a second step, we have estimated the uncertainty of the 
fitted parameters. There are standard statistical procedures 
to retrieve an uncertainty for parameters determined by a 
least-square fit. These uncertainties relate to the quality of 
that particular fit, but do not provide a meaningful estimate 
for the error bars to put on physical quantities in our case. 
Indeed, as we perform fits on average spectra, the dispersion 
of the physical parameters among the individual spectra has 
been removed. In order to get a better estimate of the error 
bars, we have fitted all individual spectra, using as seed 
values the parameters obtained for the corresponding aver- 
age spectrum. These latter fits were not "educated": each 
model parameter was left free to vary. As an illustration, 
Fig. 6 displays all individual spectra with an experimental 
gap of 44 meV, together with the fitted curves. From the 



D. Fitting procedure and accuracy of the parameters 

The 13-parameter landscape is too complex for a blind 
brute-force fitting approach. It contains plenty of local min- 
ima, providing equally satisfactory fits "by the eye", but 
parameters — especially band parameters-which are signifi- 
cantly different. The reason is that the theoretical spectrum 
does not depend critically on each hopping amplitude t t 
separately, but rather on global properties of the nonin- 
teracting DOS, such as the position and asymmetry of the 
VHS, that are not in one-to-one correspondence with the 
set of t, values. It is indeed possible to find different sets of 
hopping amplitudes, which give rise to noninteracting DOS 
curves that are hardly distinguishable in a limited energy 
window. 66 Hence the probability is close to one, that a least- 
square fit, starting with random values of the parameters, 
ends in some unphysical local minimum. It is even very 
likely that the absolute minimum, corresponding to the best 
fit in a mathematical sense, would provide for some param- 
eters, values which must be excluded on physical grounds. 
Therefore, a pinch of physical common sense must be added 
to the procedure, and the fits can not be entirely automatic. 
Such "educated fits" are rather (human) time consuming, 
and can hardly be undertaken for a large number of spectra. 



i 1 1 1 r 




-100 100 

Sample bias (mV) 



FIG. 6. Series of local experimental spectra with a peak-to-peak 
gap A p = 44 meV (yellow circles), and their fit to the microscopic 
model (red lines). Similar results are obtained for all series in 
Fig. 3. 
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distributions of fitted values in such series, we deduce an 
average and a standard deviation for each fitted parameter, 
which can be compared with the values obtained by fitting 
the average spectra. 



IV RESULTS 

Figure 7 presents the results of our fits to the average 
spectra, and Table I lists the complete set of fitted param- 
eters, as well as the estimated uncertainties. Considering 
the simplicity of the model — with regard to the high chemi- 
cal and electronic complexity of the material — the quality 
of the fits in Fig. 7 is remarkable. The model reproduces 
the relative values of the conductance at zero bias, on the 
coherence peaks, and on the background, as they vary with 
increasing A p . It also captures the electron-hole asymmetry 
of the coherence peaks, and the decrease of this asymmetry 
with increasing A p . The electron-hole asymmetry of the dip 
can be followed as a function of A„, even when the coher- 
ence peaks have become almost symmetric, at A p = 54 meV 
The theoretical spectra are systematically too "V"-shaped 
near zero bias, as compared to the experimental ones that 
look more rounded, particularly for smaller A p values. This 
might be an indication of a small deviation from pure d x 2_ y 2- 
wave symmetry close to the gap nodes. Although we have 
restricted the fits to the low-energy region, it is also possible 
to obtain good fits in the full energy range of Fig. 3, at 
the price of a slight deterioration of the fit quality at low 
energies, especially for the largest gaps. The progressive 
inadequacy of the model for increasingly extended energy 
ranges can receive different explanations. One possibility is 
the presence of additional components in the experimental 
spectrum, not considered in the model, like the contribution 
of the inner-layer band and/or a correlation-induced ex- 
cess of spectral weight for occupied states, which develops 
as A p increases. 40 Another possibility is the existence of 
other scattering mechanisms at high energy, in particular 
the scattering on the continuum of spin fluctuations. 37 

Before discussing the parameters reported in Table I, we 
emphasize three implications of Fig. 7. The excellent quality 
of the fits substantiates our basic assumption, that scanning 
tunneling spectroscopy measures the local electron DOS: 
the raw STS data can be directly compared with a model 
calculation based on physically plausible parameters. The 
presence of the Van Hove singularity in the noninteracting 
DOS is crucial to reproduce the various asymmetries of the 
spectra: the noninteracting dispersion is the only source of 
electron-hole symmetry breaking in this model. Finally, a 
coupling to the spin resonance can quantitatively explain 
the redistribution of spectral weight around the dip energy. 
The tight localization of the resonance around (n, tz), as 
opposed to optical phonons that span the whole Brillouin 
zone, plays a key role in producing the correct line-shape of 
the dip. 

Inspecting Table I, one first notices that the values ob- 
tained by fitting the average spectra coincide, within the 
uncertainty, with the average of the values obtained by 



fitting individual spectra. One also notices the large un- 
certainties, in some cases up to 100%, associated with the 
band parameters t h illustrating the weak sensitivity of the 
theoretical spectrum to each individual t ; (see Sec. HID). 
The nearest-neighbor hopping t 1 varies between —140 and 
—256 meV: these values are plausible, and fall within the 
range of published values for Bi-based cuprates. 37,58 ' 69 Phys- 
ically, one does not expect large local variations of t x , such 
as those seen in Table I. It should be noted, however, that 
these variations are comparable with the typical uncertainty 
of t 1 . Furthermore, the set of fitted hopping parameters 
should not be regarded as an accurate determination of the 
microscopic hopping amplitudes, but rather as an effective 
representation of the low-energy dispersion. We will see 
below that the properties of the dispersion which have rele- 
vance close to the Fermi surface, such as the VHS energy and 
the Fermi velocity have comparatively smaller uncertainties. 
The hopping amplitude t 3 to the second neighbor along the 
Cu-Cu axes is negative, like t 1 , as expected from symmetry 
considerations. Likewise, the diagonal hoppings t 2 and t 5 
are both positive, or the latter is almost zero (we consider 
the small negative values of order 1 meV as insignificant). 

The phenomenological scattering rate T remains small, 
with a tendency to increase for larger gaps. This variation 
allows us to capture part of the overall broadening of the 
spectra observed in the raw data with increasing A p . An- 
other part of this broadening, for A p ^ 44 meV, is provided 
by a steep increase of the spin resonance width r s . Values of 
~ 15 meV for T s seem somewhat too large, when compared 
with available values for Bi-2212 and Bi-2223, 13 ' 63 but we 
note that their uncertainty is also quite large. The uncer- 
tainty of T is comparatively smaller, because this parameter 
provides the only broadening mechanism at energies below 
f2 s , and it is therefore well constrained by the line-shape 
around zero bias. We find that A is 20 to 30% larger than 
A p . As discussed in Sec. Ill C, a reduction of the bare gap 
occurs when Aq exceeds a critical value, and the pairing 
induced by the spin resonance is over-compensated by the 
downward renormalization of the energy levels. In our 
case, Aq shows a tendency to increase with increasing A p , 
and correspondingly, the difference between A and A p 
also increases. The fitted fi 5 are anti-correlated with A p , 
as discussed in more detail below, and consistently with 
previous studies. 31 ' 32 Lastly, the dimensionless coupling a 
exhibits no clear trend, but the product (at-^ 2 , which con- 
trols the strength of the scattering rate [see, e.g., Eq. (14)], 
increases steadily with increasing A p . We discuss the cou- 
pling strength further below. 



V DISCUSSION 

A. Electron-hole asymmetry and Van Hove singularity 

The tunneling spectra exhibit systematic asymmetries be- 
tween positive and negative bias, suggestive of an intrinsic 
electron-hole asymmetry of the electron DOS. The asymme- 
tries concern the height of the coherence peaks, the strength 
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FIG. 7. Average spectra of Fig. 3 (circles) and fits with the microscopic model described in Sec. Ill (red curves). The model parameters 
are given in Table I (first column for each parameter) . 



TABLE I. Values of the model parameters obtained from fitting average and individual spectra. For each experimental gap A p and each 
parameter, the first column gives the value obtained by fitting the average spectrum, as shown in Fig. 7; the second column gives the 
average of the values obtained by fitting all individual spectra, as illustrated in Fig. 6; the number in parentheses is the corresponding 
standard deviation. 



(meV) 


/x (meV) 


ti 


Parameters of the noninteracting dispersion 
(meV) t 2 (meV) t 3 (meV) t 4 


(meV) 


t 


5 (meV) 


36 


-165 


-187 (56) 


-180 


-205 (43) 


32 43(24) 


-11 


-6(10) 


2.2 


-2.7 (8.6) 


-1.6 


-0.7 (3.7) 


38 


-125 


-164(111) 


-161 


-192 (79) 


21 31(36) 


-14 


-11(17) 


5.8 


3.0 (9.3) 


-1.4 


-3.6 (6.3) 


40 


-101 


-135(105) 


-140 


-176 (85) 


3 21(46) 


-23 


-15(21) 


11.9 


3.7(15.5) 


-3.4 


-2.3 (4.1) 


42 


-116 


-122 (51) 


-162 


-179 (46) 


18 24(26) 


-17 


-15(13) 


4.3 


-0.8 (9.2) 


-0.8 


2.9 (3.5) 


44 


-237 


-274(108) 


-206 


-247(217) 


56 52(50) 


-36 


-64(93) 


-10.3 


-13.3(39.9) 


27.9 


42.3 (42.0) 


46 


-306 


-243 (73) 


-179 


-169 (66) 


59 29(43) 


-79 


-93 (49) 


-22.6 


-19.4(12.6) 


56.6 


55.1(16.6) 


48 


-340 


-329(108) 


-206 


-206(109) 


60 36(52) 


-97 


-129(55) 


-26.3 


-20.6 (9.3) 


67.0 


76.1(24.9) 


54 


-305 


-296(146) 


-256 


-239 (87) 


58 18(44) 


-88 


-134(73) 


-35.4 


-27.1(15.9) 


62.9 


69.7(25.6) 



Scattering rate, BCS gap, and (n, n) resonance 





r (meV) 


A (meV) 


Q s (meV) 


T s (meV) 


Aq (a" 1 ) 


a 


36 


2.0 2.3(0.7) 


42.7 42.3(1.0) 


36.8 36.8(0.9) 


1.6 1.8(1.1) 


1.39 1.34(0.11) 


0.73 0.67(0.14) 


38 


2.1 2.5(0.6) 


45.0 46.0(3.7) 


34.5 34.6(2.5) 


0.5 1.0(1.1) 


1.46 1.45(0.34) 


0.79 0.76(0.19) 


40 


1.0 1.9(0.9) 


47.9 46.7(3.4) 


34.4 35.7(2.3) 


1.7 2.2(2.1) 


1.54 1.41(0.38) 


0.97 0.88(0.29) 


42 


1.0 2.2(0.9) 


52.5 49.1(4.9) 


29.9 33.3(3.2) 


1.7 2.3(2.1) 


1.65 1.38(0.49) 


0.81 0.78(0.18) 


44 


3.6 3.9(2.4) 


48.6 47.7(7.3) 


33.7 33.5(4.5) 


4.2 4.1(5.4) 


1.15 1.20(0.60) 


0.70 0.72(0.27) 


46 


6.6 6.3(1.7) 


59.7 64.0(6.4) 


21.6 14.5(9.9) 


18.6 18.3(7.0) 


1.68 2.04(0.69) 


1.16 1.65(0.32) 


48 


7.6 7.1(1.8) 


60.4 61.7(5.9) 


23.1 18.3(9.1) 


15.1 14.1(7.2) 


1.80 1.92(0.62) 


0.97 1.31(0.38) 


54 


8.5 8.5(1.7) 


70.7 71.8(6.3) 


19.2 14.3(6.8) 


13.7 14.0(8.0) 


2.17 2.42(0.58) 


0.79 1.13(0.35) 



of the dip feature, and the conductance background. The 
three effects are present, with varying magnitudes, in all 
spectra, as highlighted in Fig. 8. In the model, any asymme- 
try can be traced back to the noninteracting band structure. 
In particular, there is a simple relationship between the 



width and asymmetry of the dip, and the energy £ M of the 
VHS in the noninteracting band, 31 as is clear from Fig. 4(a). 
A wider and/or stronger dip at negative bias implies £ M < 0. 
Furthermore, one generally finds that the negative-energy 
coherence peak is taller when £ M < 0, although this rule 
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FIG. 8. Average experimental spectra of Fig. 3, with the energies 
measured from the maximum of the negative-bias coherence peak. 
The dashed lines show the negative-bias spectrum, mirrored at 
positive energies to highlight the electron-hole asymmetry. The 
circles indicate the dip minimum, the triangles show the inflection 
point (peak in d 2 I/dV 2 ), and the straight lines are guides to the 
eye. 



is not absolute [see Fig. 4(a)]. In this respect, the series 
of experimental spectra in Fig. 8 present two trends, which 
seem to have conflicting implications. On the one hand, the 
coherence peaks become more symmetric with increasing 
A p , suggesting that £ M is negative for small gaps, and ap- 
proaches zero for larger gaps. On the other hand, there is a 
tendency for the dip to become wider and more asymmetric 
with increasing A p , indicating that the VHS moves to lower 
energies with increasing A p . 

The latter evolution is consistent with the interpretation 
that spectra with larger gaps correspond to more under- 
doped regions with higher electron densities. It is also sup- 
ported by the fits, which yield larger negative values of £ M 
for larger gaps [see Fig. 9(a)]. Interestingly, by extrapolat- 
ing our results, one expects a VHS at positive energy for gap 
values lower than 24 meV, i.e., on the strongly overdoped 
side. We are not aware of any systematic investigation of 
the VHS by angle-resolved photoemission (ARPES) in Bi- 
2223. One point of comparison is given by the analysis of 
Ref. 70, which leads to a dispersion approaching —25 meV 
at (7i, 0), the value that we obtain for gaps close to 45 meV 
However, this analysis neglects renormalization effects, and 
may therefore underestimate £ M . 



B. Spin-resonance energy 

Figure 9(b) shows the fitted values of Q. s . For the smallest 
gaps, corresponding to spectra with tall and asymmetric 
coherence peaks, we find values between 30 and 40 meV, 
and a good correspondence between the Q s obtained from 
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A p (meV) A p (meV) 

FIG. 9. Physical parameters, as determined by the fits. The 
full symbols with error bars show the average parameters and 
their standard deviation, obtained by fitting all individual spectra. 
The empty symbols show the parameters obtained by fitting the 
average spectra. The crosses in (b) represent the peak-to-dip 
energy difference E A — A p of the experimental average spectra, 
corresponding to the circles in Fig. 8. (c) Renormalization factor, 
Eq. (19), at the M point (tt.O). (d) Nodal velocity. 



the fits, and the peak-to-dip energy difference, the latter 
being a few meV larger than the former. Both magnitudes 
follow the same decreasing trend with increasing A p . We 
stress that the spectral feature corresponding to a peak in the 
second derivative spectrum d 2 I /dV 2 occurs at a significantly 
higher energy of 56 meV (triangles in Fig. 8), similar to 
the values observed in Bi-2212 (Ref. 71). Therefore, the 
peak in d 2 I/dV 2 does not provide a good estimate of fi s . 
This is a hallmark of the d-wave gap symmetry, and occurs 
irrespective of the nature of the collective mode. 36 ' 39 For 
gaps larger than 45 meV, the fitted fi s drop abruptly to values 
close to 20 meV, in contrast to the peak-to-dip energies in 
the experimental spectra, that stay close to 35 meV We 
presume that part of this drop is driven by changes in the 
tunneling spectra, which go beyond the scope of the model: 
the rapid drop of fi s is accompanied by other changes, such 
as a sudden raise of T s and t 5 (see Table I), suggesting that 
the fit has been pushed to a different region of the parameter 
space. 

Inspecting the data, one sees that an obvious change 
between the 44 and 46 meV spectra is a falloff of the peak 
to background ratio. This modification cannot be accounted 
for by simply increasing F and/or F s , since this would also 
reduce the dip to background ratio, while in the spectra the 
latter remains unchanged. The compromise found by the 
fitting routine is to lower Q s , and thus to reduce the peaks' 
height without affecting the dip strength [see Fig. 4(b)], 
and to tune the position of the dip minimum by adjusting 
other parameters [in particular, increasing Aq pushes up 
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the ratio A /A p , and increases the peak-to-dip separation 
as seen in Fig. 4(c)]. Another trend is that the coherence 
peaks become symmetric for large gaps, and this also drives 
Q s downward in the fits. 

In the absence of neutron scattering data for fi 5 in Bi- 
2223, and considering the large uncertainty in the fitted 
values for A p larger than 45 meV, it is difficult to ascertain 
whether the actual decrease of the spin-resonance energy 
for large gaps is as sudden as suggested by our fits, or 
rather more continuous. Nevertheless, the range of values 
which we find corresponds well to the doping evolution 
observed in Y-123, where the resonance energy decreases 
from ~ 40 meV at optimal doping to 25 meV in under- 
doped samples. 72 ' 73 In Bi-2212, the resonance is found at 
42 meV at optimal doping, and its energy decreases down 
to 34 meV with strong overdoping, 13 but no data has been 
so far reported in the underdoped region. 

C. Coupling strength and renormalization factors 

The values of the bare coupling strength obtained from 
our fits, i.e., g in Eq. (5), are similar to those used in Ref. 37. 
However, as discussed in Sec. Ill, comparing the bare cou- 
plings obtained by fitting different spectra can be misleading, 
because the strength of the strong-coupling effects depends 
also on other model parameters, particularly on the momen- 
tum spread Aq of the resonance. A physically meaningful 
measure of the coupling strength is provided by the renor- 
malization factor, which we define as 

^k = ~~ [£n(k ) 60 + E 2 2(*,«<0] a>=0 . (19) 

0.CO z 

The Brillouin-zone average of this quantity is given approx- 
imately by Eq. (15). The evolution of A fc at the M point 
with varying A p is shown in Fig. 9(c): A M increases sharply 
at A p = 45 meV Like for fi s , the A M error bars are larger 
for the spectra with larger gaps. The renormalization A fc is 
anisotropic in momentum space, as shown in Fig. 10. The 
values of A fc are smaller in the nodal region than close to 
the anti-nodal points. The figure also shows that, for the 
smallest gaps, the Fermi surface average of A fc (dots) is well 
approximated by the full Brillouin-zone average (horizontal 
bars) given by Eq. (15). For larger gaps, the Fermi surface 
moves away from the M point, where the renormalization 
is strongest and, as a result, A is larger than the average 
Fermi-surface renormalization. We note that the expression 
(17), valid for Aq = 0, overestimates A M by a factor ~ 2.5 
for small gaps, where Aq is smallest, and by a factor ~ 3.5 
for the larger gaps with larger Aq values. The formula (17) 
should therefore be used with caution, unless Aq is very 
small. The fast increase of A M above A p = 45 meV is ac- 
companied by an increase of anisotropy: the difference in 
renormalization between anti-nodal and nodal regions is 
20-30% for the smaller gaps, and increases to 30-50% for 
the larger gaps. Note that the X k shown in Fig. 10(a) are 
significantly larger, and more anisotropic, than the values 
found for Bi-2212 using phonon models. 74 This suggests that 
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FIG. 10. Renormalization factor A fc on the renormalized Fermi 
surface, (a) Superconducting state: the parameters are those 
corresponding to the average spectra of Fig. 7 with the corre- 
sponding value of A p written on each curve. The dots indicate 
the Fermi-surface average of X k , and the horizontal bars show 
the Brillouin-zone average calculated with Eq. (15) (off-scale for 
A p = 46 meV). Inset: Fermi surface for A p = 42 meV, and defini- 
tion of the Fermi-surface angle cp. (b) Normal state: the solid lines 
show the value of k k obtained by setting A = and T = 110 K, 
and keeping the other parameters unchanged. The dashed lines 
correspond to A = 0, T = 110 K, and r s = 25 meV 



fitting our STM data to phonon-based models would lead to 
unrealistically large electron-phonon matrix elements. 

The experimental determination of A by ARPES and op- 
tical conductivity are often performed in the normal state. 
In view of a comparison with those data, we note that the 
calculated A's increase significantly on going from the su- 
perconducting to the normal state. This can be understood 
by looking at Eq. (15), which shows that A increases as 
the gap in JV (e) closes. Figure 10(b) displays the A values 
computed at the Bi-2223 critical temperature T c = 110 K, 
with A set to zero, while keeping the other parameters un- 
changed. The renormalization increases by a factor 1.2-2.4 
with respect to its value in the superconducting state, and 
the anisotropy also raises up to 30-60%. Experimentally, 
the transition to the normal state also affects the spin reso- 
nance, which broadens in energy 75 on warming across T c . 
This can be modeled by increasing the value of T s , while 
keeping the other parameters fixed. In the absence of ex- 
perimental input on the normal-state value of F s in Bi-2223, 
we show in Fig. 10(b), as an illustration, the normal-state 
renormalization calculated with T s = 25 meV The values of 
the renormalization are reduced by the increase of r s , and 
are similar to the values calculated in the superconducting 
state if T s = 25 meV 

Our renormalization factors compare well with experi- 
mental values reported in the literature. In underdoped 
and overdoped Bi-2212, the normal-state renormalization 
near the anti-nodal point, determined by photoemission, 
was found to be 1.5. 26 This is consistent with our finding 
A M fcf 2, considering that the more pronounced dip feature 
in the STM spectra of Bi-2223 with respect to Bi-2212 indeed 
points to a larger coupling in the three-layer compound. In 
the normal state, but at the nodal point, a renormalization 



decreasing from 0.8 to 0.55 as a function of increasing hole 
doping was measured at T = 120 K in Bi-2212. 76 Our nodal 
values for Bi-2223 in the normal state, and for T 5 = 25 meV, 
decrease from 0.8 to 0.5 with decreasing gap size, which 
draws a similar trend as observed in Bi-2212. Below T c , 
nodal values of A between 0.7 and 0.9 for Bi-2212 are re- 
ported in Ref. 76, very close to our superconducting-state 
results of Fig. 10(a). In photoemission experiments, an 
accurate determination of the anti-nodal renormalization 
at low temperature is difficult, due to the presence of the 
superconducting gap. A detailed study reported a value of 2 
for optimally-doped Bi-2212, 77 while our values for Bi-2223 
vary between 1 and 1.6. To our knowledge, no ARPES value 
of X has been explicitly reported for Bi-2223. Nevertheless, 
from a recent study, 78 we may estimate a nodal A of 0.6 at 
10 K for an optimally-doped sample, which lies at the lower 
edge of our values in Fig. 10(a). Finally, the average renor- 
malization factor in the normal state of optimally-doped 
Bi-2223 has been estimated by fitting optical data 35 to a 
model involving a spectrum of bosonic excitations, and lead 
to the values 2.18 and 1.75, depending whether the full 
bosonic spectrum, or only its low-energy part, is taken into 
account, respectively. Both values fall within our range of 
normal-state renormalizations in Fig. 10(b). 

Figure 9(d) shows the nodal velocity, calculated using the 
parameters obtained from the fits, and assuming a Bi-2223 
in-plane lattice constant of 3.825 A. The nodal velocity can 
be directly measured by ARPES, unlike the renormalization 
X k , which requires assumptions about the noninteracting 
dispersion. Ref. 78 reports nodal velocities between 1.5 
and 1.7 eV A for the outer-layer band of optimally-doped 
Bi-2223, i.e., 2.3-2.6 x 10 7 cm/s. These values agree well 
with our results for A p > 44 meV In the case of Bi-2212, 
values ranging between 1 and 2.5 x 10 7 cm/s as a function 
of doping have been reported. 79 We conclude that the varia- 
tions seen in Fig. 9(d) are consistent with variations of the 
nodal velocity associated with changes in the local doping 
level. 



D. Spectral function and simulated ARPES intensity 

The parameters determined by fitting STM spectra allow 
us to calculate the momentum-resolved spectral function, 
and to make predictions for the measured ARPES inten- 
sity The generic properties of the spectral function in the 
present model have been extensively discussed. 20 ' 37 The 
main characteristics are summarized in Fig. 11. At energies 
below the onset of scattering, |w| < f2 s , the dispersion is 
renormalized downwards, but not broadened. This can be 
seen most clearly along the nodal direction (0, 0)-{n, ri) in 
Fig. 11(a). At these low energies, the quasiparticles form 
banana-shaped regions around the nodal points [Fig. 11(b)] . 
Increasing energy from the onset of scattering at fi s , the scat- 
tering rate increases, and reaches its maximum at A + Q s 
[see Eq. (14)]. Since A p <A G + fi s , the anti-nodal excita- 
tions at \co\ = A p remain rather sharp [Fig. 11(c)]. When 
the scattering rate is maximum at \co\ = A + fi s , in con- 
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FIG. 11. Spectral function calculated using the parameters cor- 
responding to the average spectrum with A p = 42 meV in Fig. 7. 
(a) Along high-symmetry lines in the Brillouin zone. The color 
scale shows the variation of the spectral function from zero (white) 
to its maximum (yellow). The orange line is the noninteracting 
dispersion. The red line shows the BCS dispersion; the width of 
the line is proportional to the spectral weight, (b), (c), and (d) 
Cuts at three characteristic energies. The color scale is the same in 
all graphs. 

trast, the excitations are very broad [Fig. 11(d)]. At this 
energy, the real part of the self-energy changes sign, and the 
quasiparticle dispersion correspondingly jumps from below 
to above the noninteracting dispersion. This abrupt change 
in dispersion renormalization occurs at the same energy in 
the whole Brillouin zone, both for occupied and empty states, 
and leads to a suppression of states which is responsible 
for the dip in the DOS. Near the anti-node, the d-wave gap 
induces additional structures: the minimum-gap locus is 
close to — but not exactly at — the position of the noninter- 
acting Fermi surface along the (71, 0)-(7i, n) line, due to 
Fermi-surface renormalization. One also notices the reduc- 
tion of the gap A p , as compared to the BCS gap A , due 
to the competition between pairing and renormalization, 
as discussed in Sees. Ill C and IV. In the low-energy region, 
the weakly-dispersing quasiparticle branch near (n;, 0) has 
lower energy than the corresponding BCS branch, while 
above the dip energy, the quasiparticle energy is higher than 
the noninteracting and BCS dispersions, like in the nodal 
region. 

In order to quantitatively compare our results with real 
ARPES data, we have used our model to simulate the upper 
panels of Fig. 2 in Ref. 80. With the aim of reproducing 
the experimental conditions of this work, we set the tem- 
perature to 10 K, and filter our spectral function with a 
Gaussian representing an energy resolution of 18 meV, and 
a momentum resolution of 0.05/a. 80 The result is displayed 
in Fig. 12. Close to the nodal direction (left and middle 
panels), the agreement is good. The model — which only 
includes the spin resonance — is too crude to completely cap- 
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FIG. 12. Simulated ARPES intensity for optimally-doped Bi-2223 at temperature T = 10 K. The spectral function of Fig. 1 1 was multiplied 
by the Fermi function, and filtered to mimic a momentum resolution of 0.05/a and an energy resolution of 18 meV The color scale and 
the three momentum cuts correspond approximately to those of Fig. 2 in Ref. 80. The black and red lines show the noninteracting and 
quasiparticle dispersions, respectively. The arrows indicate the "kink", where the quasiparticle dispersion deviates from its low-energy 
linear behavior. 



ture the measured dispersion in the region of the dip: in this 
energy range, the simulated dispersion jumps from below to 
above the noninteracting dispersion, while the experiment 
interpolates smoothly across the jump. We attribute this dis- 
crepancy to additional scattering mechanisms, not included 
in our model, in particular those involving the continuum 
of spin fluctuations. 37 

The energy of the "kink", around — 70 meV, is well repro- 
duced by the calculation. We note, however, that the feature 
corresponding to the energy A + fi s is not the kink, but the 
midpoint of the jump, where the quasiparticle dispersion 
crosses the noninteracting dispersion. It is also worth men- 
tioning that this energy is not dispersing, and is given by 
A +fi s , in spite of the fact that the gap vanishes at the node, 
as emphasized in Fig. 11(a). The experimental determina- 
tion of this energy scale by ARPES requires an assumption 
for the shape of the noninteracting dispersion. Our results 
call into question the assumption made in Ref. 80, namely 
that the quasiparticle and noninteracting dispersions meet 
near —200 meV This assumption has direct implications on 
the value of the self-energy derived from ARPES data, in 
particular, on the energy at which the real part vanishes. 
The value near —200 meV reported in Ref. 78 is obviously 
a direct consequence of the assumption made for the non- 
interacting dispersion. Using the dispersion resulting from 
our fits, one would be lead to the conclusion that the real 
part of the self-energy derived from ARPES vanishes at the 
energy — (A + fi s ). In the anti-nodal region (right panel 
of Fig. 12), there are differences between our results and 
the ARPES data: the low-energy part, below the kink, is 
too dispersive in the model, while the high-energy part is 
not dispersive enough. This could be partly due to different 
Fermi surfaces in the experiment and in the model, which 
imply that the segments of the dispersion considered in the 
experiment and in the model are not exactly identical. It 
could also be the consequence of high-energy scattering 
processes, neglected in the model. Despite these differences, 
the kink energy and the minimum in the spectral intensity, 
between the low-energy and high-energy parts, are very 
similar in the simulation and in the experiment. The overall 



agreement between Fig. 12 and Fig. 2 of Ref. 80 support 
the claim that STM tunneling spectra, although they come 
from a momentum-integrating probe, do contain the nec- 
essary information needed to reconstruct the low-energy 
momentum-resolved spectral function. 



VI. SUMMARY AND CONCLUSIONS 

We have performed an analysis of STM spectra measured 
on optimally-doped Bi-2223, by means of a microscopic 
model which takes into account the Van Hove singularity 
within a one-band description, a BCS gap with pure d-wave 
symmetry, and a coupling to the (71, n) spin resonance. This 
model can reproduce the experimental spectra remarkably 
well (Fig. 7), and provides values for the relevant physical 
parameters that are sound, and consistent with values ob- 
tained by other experimental probes. The inhomogeneity 
of the electronic properties observed on the sample surface 
allowed us to study variations in parameters as a function 
of the spectral gap A p . Assuming a one-to-one correspon- 
dence between A p and the hole doping level, changes in A p 
may be interpreted as variations of local doping. The main 
trends are that the Van Hove singularity moves to lower 
energies, and the energy of the spin resonance decreases, 
with increasing A p (Fig. 9). The former supports the in- 
terpretation that larger gaps correspond to lower dopings, 
and the latter supports the claim that the dip is caused by 
the spin resonance, rather than by optical phonons. The 
strength of the coupling to the spin resonance, measured by 
the dispersion renormalization, increases steadily for gaps 
larger than 42 meV, namely towards the underdoped region 
of the phase diagram. 

The presence of a Van Hove singularity, breaking the 
electron-hole symmetry of the electronic spectrum, is un- 
mistakable in the raw data (Fig. 8). Nevertheless, the pre- 
cise determination of the hopping amplitudes t i remains a 
challenge, as illustrated by the large uncertainties attached 
to these parameters in Table I. The reason is that relatively 
large variations of the t, 's can collaborate to induce marginal 
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changes in the electron DOS, which largely govern the STM 
spectrum. Our confidence in the fitted tight-binding disper- 
sion stems from the ability of the whole model to reproduce 
momentum-resolved ARPES data with good accuracy. In 
turn, our determination of the band-structure may provide 
indications on how to extract the self-energy from ARPES 
measurements. 



ACKNOWLEDGMENTS 

We acknowledge useful discussions with M. Eschrig. This 
work was supported by the Swiss National Science Founda- 
tion through Division II and MaNER 



Appendix A: The function B{co,E) 

Equations (5) to (11) imply that the function B(co,E) 
entering Eq. (13) reads 



B{co,E) = A 2 de l de 2 [L T {e l -tl s )-L T {£ 1 + n s )} 

J — 00 

b(gx)+/(-g 2 ) .... 

x L f [e 2 - £) — — , (Al) 

co — e ± — e 2 + iCr 

where b and / are the Bose and Fermi functions, respec- 
tively. Using the identity JdxL r (x)/(z — x) = l/[z + 
irsign(Imz)], we obtain 
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The symbol in braces means that the result obtained from 
the first two terms on the r.h.s. must be anti-symmetrized 
with respect to fi s . The remaining integrals can be evaluated 
by closing the integration contour in the complex plane. The 
integrands can we rewritten in terms of products of simple 



poles by means of the identities 
1 f 1 

L r ( Z ) : 
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(A3) 
(A4) 

(A5) 



Only one half of the poles of the Bose or Fermi function are 
enclosed in the contour and give a contribution. Therefore, 
the integrals involve semi-infinite sums over Matsubara fre- 
quencies, which can be converted into the digamma function 
ip, using the relation 



The final result is 
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The function B simplifies considerably in the case of a sharp 
resonance (T s = + ) and for sharp Bogoliubov quasiparticles 
(f = + ), as well as zero temperature. In this case, we have 



A 2 



co-E-n s sign(£) + i0+ ' 



(A8) 



as can be readily deduced from Eq. (Al), by replacing the 
Lorentzians by delta functions. 
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